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Abstract — A very strong and persistent phytoplankton bloom was observed by ocean color 
satellites during September - December 2003 along the northern Patagonian shelf. The 2003 bloom 
had the highest extent and chlorophyll a (Chl-a) concentrations of the entire Sea-Viewing Wide 
Field-of-view Sensor (SeaWiFS) period (1997 to present). SeaWiFS-derived Chl-a exceeded 20 mg 
m‘ 3 in November at the bloom center. The bloom was most extensive in December when it spanned 
more than 300 km across the shelf and nearly 900 km north - south (35° S to 43° S). The northward 
reach and the deep penetration on the shelf of the 2003 bloom were quite anomalous when compared 
with other years, which showed the bloom more confined to the Patagonian shelf break (PSB). The 
PSB bloom is a conspicuous austral spring-summer feature detected by ocean color satellites and its 
timing can be explained using the Sverdrup critical depth theory. Based on high-resolution numerical 
simulations, in situ and remote sensing data, we provide some suggestions for the probable 
mechanisms responsible for that large interannual change of biomass as seen by ocean color 
satellites. Potential sources of macro and micro (e.g., Fe) nutrients that sustain the high 
phytoplankton productivity of the Patagonian shelf waters are identified, and the most likely physical 
processes that maintain the nutrient balance in the region are discussed. Seasonal and interannual 
changes in vertical mixing, shelf intrusions of the Malvinas Current, and the position of the Brazil- 
Malvinas confluence zone, are the most likely mechanisms driving the observed changes in 
phytoplankton production. Previous studies advocate that the South Atlantic Ocean is the recipient of 
dissolved Fe from various sources, including, atmospheric deposition from desert dust and volcanic 
ash, and coastal submarine groundwater discharge from lagoon systems, but the scope of our 
analysis and data limitations precluded us from establishing a connection between these sources of 
Fe and the Patagonian shelf bloom. This paper highlights the complexity of the biogeochemical 
processes that control the shelf biomass and production of the region and provides suggestions for 
further studies to improve our understanding of the physical-biological interactions responsible for 
nutrient balance and ecosystem response. 
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1. Introduction 


Along the Patagonian coast off Argentina, the shelf-break front formed between shelf waters 
and the Malvinas (Falkland) Current (MC) shows a conspicuous band of high phytoplankton 
biomass throughout spring and summer detected by ocean color sensors. This region is the feeding 
and spawning ground of several commercial species of fish and squid [Acha et al., 2004] and is 
thought to play an important role in CO 2 sequestration by the ocean [Bianchi et al, 2005], The 
physical and biological characteristics of the Patagonian Shelf Break (PSB) front have been 
previously investigated [Acha et al., 2004; Rivas, 2006; Rivas et al, 2006; Saraceno et al, 2005; 
Signorini et al, 2006], The region bordering the PSB presents elevated chlorophyll a (Chl-a) from 
spring through autumn presumably supported by the nutrient-rich subsurface waters of the MC. This 
nutrient supply is probably a result of upwelling processes caused by the interaction between the MC 
and the continental slope [ Carreto et al, 1995; Matano and Palma, 2008], Near 39°S there is an 
energetic southward and eastward flow resulting from the collision of the southward Brazil Current 
(BC) with the MC, namely the Brazil/Malvinas Confluence (BMC). The MC transports cold and 
relatively fresh sub-Antarctic waters equatorward, and the BMC generates one of the most energetic 
regions of the world ocean [Piola and Matano, 2001], The BMC extends eastward forming the 
South Atlantic Current (SAC). The SAC is an eastward conduit of biogenic material originating on 
the outer shelf and slope through meanders and eddies, as can be seen from satellite images (Figure 
1 ). 

It has been suggested that low productivity in nutrient-rich surface waters of the Southern 
Ocean may result from a deficiency of Fe [Martin et al, 1990], The MC is a branch of the Antarctic 
Circumpolar Current (ACC) and thus inherently low in iron concentration. There is in situ evidence 
that the Malvinas upwelling region is a large source of macro nutrients, especially nitrate [Carreto et 
al, 1995; Garcia et al, 2007], but the origin of Fe to maintain the intense PSB bloom remains 
elusive due to the sparse in situ Fe measurements. Bowie et al [2002] made a few measurements of 
Fe in the region and reported a maximum of 1.7 nM in the front of the sub-tropical convergence 
zone towards Uruguay at 35°S. Further south, they found that the total dissolved Fe tended towards 
baseline values (~0.3 to 1.0 nM), with an increase (2.7 nM) observed in the shelf waters off the East 
Malvinas Island, but there are very few samples south of 35°S. Two sources of Fe for the surface 
waters of the Atlantic Southern Ocean (SO) have been proposed [Cassar et al, 2007; Erickson et al, 
2003; Gaiero et al, 2003; Meskhidze et al, 2007]: (1) oceanic input from upwelling and lateral 
flows from continental margins; and (2) atmospheric input from deposition of mineral dust 
emanating from the arid regions of South America and ash clouds from volcanic eruptions. 
Regardless of its origin, changes in Fe concentration can significantly alter ocean productivity. With 
increasing dissolved Fe concentrations, the growth rates of open-ocean SO diatoms can increase 
several fold [Timmermans et al, 2004]. 

In this paper we provide evidence of large interannual variability of ocean productivity along 
the Patagonian shelf, with a very large and anomalously intense bloom occurring in spring-summer 
2003. Based on high-resolution numerical simulations and remote sensing data, we provide some 
suggestions for the probable mechanisms responsible for that large interannual change of biomass as 
seen by ocean color satellites. Moreover, we focus on the 2003 anomalous event, as it clearly stands 
out from the region’s pattern observed along more the 10 years of satellite ocean color data. 
However, the physical processes that maintain the overall balance of nutrients and resulting 
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biogeochemical interactions are not fully understood and further field and numerical studies are 
needed to better describe the observed phytoplankton variability. 
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Figure 1 . Monthly composites of SeaWiFS Chl-a for December 2003 and 2004. The corresponding 
HYCOM 30-meter currents are superimposed. The white line indicates the location of the cross-shelf 
transect used in the analysis. 

2. Data Sources and Methodology 

2.1 Satellite Data 

SeaWiFS Chl-a (production algorithm) and Photosynthetic Available Radiation (PAR) 9-km 
monthly data, and 4km Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua sea-surface 
temperature (SST) were obtained from the NASA Ocean Color distribution archive 
(ftp://oceans.gsfc.nasa.gov/') . A calcite, or particulate inorganic carbon (PIC), algorithm [. Batch et al., 
2005] was used to produce monthly 4-km PIC composites using SeaWiFS Global Area Coverage 
(GAC) data. 

Monthly composites, cross-sections and time series extracts were processed and analyzed. 
Monthly 4-km Pathfinder Advanced Very High Resolution Radiometer (AVHRR) SST 
(http://podaac-www.ipl.nasa.gov/) and monthly 0.5°x0.5° Quick Scatterometer (QuikScat) winds 
(http://www.ifremer.fr/cersat/en/index.htm') were used to derive SST time series and monthly wind 
and SST anomalies concurrent with the SeaWiFS mission. We used all available concurrent data. 
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which are limited to September 1997 - December 2006 for the Pathfinder AVHRR SST and August 
1999 - September 2007 for QuikScat winds. 


2.2 HYbrid Coordinate Ocean Model or HYCOM 

HYCOM products (mixed layer thickness, mixed layer velocity components, and sub-surface 
temperature and salinity) were obtained from the HYCOM data service web site 
(http :/ /hycom .rsmas .miami . edu/dataserver/) . Computations are carried out on a mercator grid 
between 78°S and 47°N (1/12° equatorial resolution). The horizontal dimensions of the global grid 
are 4500 x 3298 grid points resulting in ~7 km spacing on average. There are 32 vertical layers. 
Bathymetry is derived from a quality controlled NRL DBDB2 dataset. Surface forcing is from Navy 
Operational Global Atmospheric Prediction System (NOGAPS) and includes wind stress, wind 
speed, heat flux (using bulk formula), and precipitation. We obtained outputs from the hindcast 
experiment which uses the Navy Coupled Ocean Data Assimilation (NCODA) system for data 
assimilation (SST, SSH, XBTs, ARGO floats and moored buoys). The hindcast products are 
available from November 2003 until current. HYCOM is a very mature model and a vast literature is 
available on its performance, validation, and applications to a variety of ocean basins [Chassignet et 
al. , 2006; Chassignet et al., 2007; Cheng et al., 2007; Halliwell, 2004; Kara and Hurlburt, 2006; 
Kelly et al, 2007 ; Winther and Johannessen, 2006]. 

2.3 Temperature-derived Nitrate 

Figure 2a shows the station locations from which the concurrent nitrate and temperature data 
originated to derive a NO 3 vs. T relationship. Only quasi-synoptic transects that sampled the Brazil 
and Malvinas Currents and their convergence region were used. A robust regression ( 1 ^= 0 . 93) of 
nitrate versus temperature (Figure 2b) was developed based on 976 observations from 74 stations 
available from the NOAA NODC World Ocean Database (http://www.nodc.noaa.gov/OC5/). The 
regression equation is 

logio (NO3) = -0.003999 T 2 + 0.009954 T + 1 .488 (1) 

Data shallower than 50 m were not used to reduce biases due to biological consumption. Most of the 
data originated from depths between 50 and 500 m. The goal was to obtain nitrate concentrations 
resulting from physical transport only. Carreto et al. [1995] provided a linear equation (r 2 =0.97) of 
nitrate versus temperature using 32 observations from a cross-shelf transect off Mar del Plata in the 
Argentinean Sea, but we were not able to apply it to our region since it yields negative retrievals for 
temperatures greater than 12°C. 

3. Results and Discussion 

3.1 Evidence of Anomalous Bloom Behavior 

A very strong and persistent phytoplankton bloom occurred during September - December 
2003 at the northern Patagonian shelf (north of 43°S). The SeaWiFS Chl-a peaked in November and 
exceeded 20 mg m° at the bloom center (~ 40.5°S and 58°W). The bloom was widest in December 
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when it spanned more than 300 km across the shelf and nearly 900 km north - south (35° S to 43° S). 
In contrast, in December 2004 there was a much weaker bloom confined to the shelf break, which is 
a more typical behavior for the region. Figure 1 contrasts the December 2003 and 2004 SeaWiFS 
Chl-a distribution clearly showing the anomalous size and strength of the 2003 bloom. 



Figure 2. (a) NODC station locations with bathymetry (30, 100, 1000, 3000, and 3000 m) and approximate location of 
Brazil and Malvinas Currents superimposed, (b) Nitrate polynomial regression (logi 0 [N0 3 ] versus temperature) based 
on 976 in situ observations below 50 m depth. The color indicates the observation depth. 


Figure 3 shows cross-shelf time series (see Figure 1 for transect location) of SeaWiFS Chl-a 
and PAR, and surface NO3 derived from Pathfinder AVHRR SST and equation (1). Cross-shelf 
average time series are superposed in white and the horizontal black line is the position of the 100 m 
isobath (shelf break). The bloom seasonal cycle is clearly shown, with peaks during the cross-over of 
elevated PAR and relatively high NO3 occurring in November-December. Highest NO3 levels (> 20 
pM) occur in the peak of winter (July-August) in the PSB when Malvinas waters are vertically 
mixed due to convective overturning and biological consumption is at a minimum. The winter- 
summer contrast in stratification will be shown later on when HYCOM cross-shelf transects are 
discussed. 

Large interannual variability of Chl-a can be seen in Figure 3. The transect crosses the widest 
part of the December 2003 bloom and reveals that the bloom was anomalous in its strength, size, and 
position on the shelf. The normal seasonal cycle of the bloom occurs during the same period but with 
much lower Chl-a and with peak values more confined to the PSB. Figure 3 also shows that the 
spring-summer bloom transitioned to more productive conditions in 2002 and returned to more 
typical conditions in 2005. Note that the strongest bloom confined to the PSB occurred during the 
spring-summer of November 2002 - January 2003. 

The PSB is also characterized by a transition of phytoplankton groups during the spring- 
summer bloom [Garcia et ai, 2007; Signorini et al., 2006], Although very few samples have been 
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collected contemporaneously with satellite data to validate ocean color algorithms for this region, a 
time series image analysis of calcite and chlorophyll algorithms in the shelf-break region indicate 
that the blooms initiate in austral spring, probably dominated by diatoms, but in mid summer 
(December, January) the community is mainly represented by calcite-producing (coccolithophores) 
phytoplankton [Sign or ini et a!., 2006], Signorini et al. [2006] hypothesize that the transition from 
diatoms to coccolithophores occurs in mid summer when the mixed layer gets much shallower 
reducing the nutrient level at surface layers but allowing an increase of light harvesting by 
phytoplankton. The coccolithophores thrive in high light and low nutrient conditions, while the 
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Figure 3. Cross-shelf time series of (a) SeaWiFS Chl-a (mg/m 3 ), (b) SST-derived surface N0 3 (pM), and 
(c) SeaWiFS PAR (mol m 2 /d). The black contours are 3, 5, 10 mg m' 3 for Chl-a and 20 and 25 pM for N0 3 . 

The horizontal black line is the location of the 100 m isobath representing the shelf break. Time series of 
cross-shelf averages are superposed in white with proper scales on the right of each tier. 

diatoms growth is limited under these conditions and a shift from diatoms to coccolithophores 
gradually occurs. Therefore, the mixed layer thickness plays a major role in this process. 
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Figure 4. Monthly composites of HYCOM mixed layer 
thickness, SeaWiFS PIC, and SeaWiFS Chl-a for December 2003 
(left column) and December 2004 (right column). The shelf break 
is represented by the 300 m and 1000 m isobaths. 


on the shelf production as well. 


3.2 SST, Chl-a, and Wind Anomalies 


To better illustrate this hypothesis, Figure 
4 shows December 2003 and 2004 
composites of HYCOM mixed layer depth 
(MLD), SeaWiFS-derived PIC 
concentration, and SeaWiFS Chl-a. The 
PIC concentration and the coccolithophore 
flag [Brown and Podesta, 1997; Brown 
and Yoder , 1994] used in the Chl-a 
algorithm are well correlated in space and 
time [Signorini et al., 2006] suggesting 
that PIC concentration is a proxy for 
coccolithophore abundance. Note that 
there is a band of low MLD following the 
PSB. This is a result of the summer 
warming above the shallow thermocline 
along the MC upwelling zone. Although 
the accuracy of the HYCOM MLD is 
unknown for this region, it will be shown 
later (see Figure 9) that the MLD along 
the cross-shelf transect during December 
2004 was on average 7 m thinner than the 
MLD during December 2003. This is a 
large difference since the largest MLD 
along the cross-shelf transect was about 
20 m. This shallower MLD apparently had 
a positive impact on the coccolithophore 
production in mid-summer 2004 judging 
from the relatively high values of PIC (>6 
mmol C m‘j along the PSB during 
December 2004 (Figure 4). Conversely, 
the PIC concentrations along the PSB 
were much lower during December 2003, 
while the Chl-a over the shelf region 
between 35°S and 45°S were much higher 
(> 10 mg m' 3 ) than in December 2004. 
Thus the MLD variability had an impact 


To further illustrate the anomalous nature of the 2003 bloom. Figure 5 and 6 show the distributions 
of Chl-a, SST, and wind anomalies for October, November, and December 2003 - 2004, 
respectively. In contrast with the scalar parameters, caution must be taken when interpreting the 
wind anomalies (speed and velocity vectors). They are somewhat difficult to interpret as the mean of 
the speed can be quite different than the mean of the vectors. For instance, variable wind direction at 
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a particular location may have a high mean speed but almost zero mean vector velocity if the wind 
direction has a symmetric distribution. Conversely, persistent and steady winds have both elevated 
speed and vector means. Because of this, the monthly mean wind speed and vectors are also 
provided in Figure 5 and 6 as a reference. Comparison between Figures 5 and 6 shows interesting 
contrast in anomaly patterns between spring-summer of 2003 and 2004. During 2003, the shelf 
region to the south of the Rio de la Plata mouth exhibits large positive anomalies of Chl-a (up to 10 
or more mg m' 3 ), negative SST anomalies, and stronger offshore and northerly winds. During 2004 
the anomalies of Chl-a and SST reverse showing lower Chl-a and warmer SSTs on the shelf region, 
especially during December. The winds were parallel to the coast (northeasterly) blowing along the 
PSB and against the direction of the MC leading to an offshore Ekman transport component. Thus, 
coastal upwelling is not a primary factor with respect to the 2003-2004 differences. In contrast, the 
offshore winds during 2003 favored the MC northeastward flow via the Ekman transport, which is 
most probably the cause for a nearly 2° northward displacement of the BMC as shown in Figure 2. 
This displacement had an impact on the Chl-a which shows a positive anomaly extending northward 
up to the latitude of the Rio de La Plata mouth near 35°S. The 2003 cold SST anomalies are in 
agreement with larger MLDs, especially during December 2003. 

3.3 Role of Light on the Timing of the Spring Bloom and the Sverdrup Critical Depth 

Sverdrup’s critical depth hypothesis [ Sverdrup , 1953] has been applied to many oceanic and 
even estuarine ecosystems with the goal of evaluating the roles of light availability and vertical 
mixing leading to phytoplankton spring blooms [. Nelson and Smith, 1991; Obata et al., 1996; Siegel 
et al, 2002; Townsend et al., 1992], Before the iron hypothesis was introduced, the prevailing idea 
of the control of primary production in the Southern Ocean was that relatively weak vertical stability 
and strong winds combined to mix the upper ocean and its resident phytoplankton assemblages over 
depths that resulted in low time-integrated irradiance for the photosynthetic cells. Following this 
reasoning, the productivity of most of the Southern Ocean is low even in summer because the depth 
of the surface mixed layer (Zm) approaches or exceeds the Sverdrup critical depth ( Zcr ). This depth 
is a function of the solar irradiance at the sea surface, the clarity of the water, and the phytoplankton 
compensation irradiance (irradiance at which photosynthesis equals respiration rates). This 
simplified formulation assumes that the mixed layer is not nutrient-limited and that there are no 
grazing or other losses. The Zcr equation is as follows 


Zcr _ I 0 
(l - exp(-AT Zcr)) K I c 


( 2 ) 


where I 0 is the incident PAR at the surface, K is the diffuse attenuation coefficient for PAR, and I c is 
the average irradiance for the water column that equals the compensation irradiance. Sverdrup 
defined the critical depth as the depth at which the vertically integrated rates of plankton 
photosynthesis and respiration are equal. Therefore, when Zm > Zcr , the vertically integrated rate of 
phytoplankton photosynthesis is less than that needed to keep pace with consumption of organic 
carbon from respiration. Under these conditions net phytoplankton biomass cannot take place. 
Conversely, when Zm « Zcr, production exceeds respiration and the phytoplankton assemblages 
flourish. With the above criteria and assumptions in mind, we applied the critical depth analysis to 
the Patagonian shelf bloom. The value of Zcr calculated from equation (2) is subject to large 
uncertainties for several reasons. 
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In the original 
Sverdrup formulation I Q is 
taken as total solar 
irradiance rather than PAR. 
There are several other 
losses involved other than 
respiration, such as sinking, 
grazing, and release of 
dissolved organic material 
(DOM). Although not all 
uncertainties can be 
reduced or eliminated, 
Nelson and Smith [1991] 
reformulated equation (2) 
by substituting terms that 
are more consistent with 
‘present-day’ optical and 
biological data while still 
preserving the Sverdrup’s 
original logic and 
mathematics. Using their 
reformulated equation, they 
investigated the control of 
Southern Ocean (including 
the ACC) productivity by 
the irradiance-mixing 
regime. The reformulated 
equation for Zcv is 
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Figure 5. Monthly anomalies of SeaWiFS Chl-a (A, B, and C), AVHRR SST (D, 
E, and F), and QuikScat wind speed and vectors (G, H, and I) for October, 
November and December 2003. QuikScat absolute monthly wind and vectors are 
shown at the bottom (J, K, and L). 


where the coefficient 3.78 
is a result of a series of 
transformations, including 
a 20% surface PAR 
reduction due to surface 
reflectivity, but is the 
equivalent of I c in units of 
mole quanta m' 2 day' 1 , 
which are the same units of 


I 0 . K is in m' 1 and Zcr in m. They also used the following empirical relationship first proposed by 
Riley [1956] for K dependence on chlorophyll for systems where most absorbance of irradiance is by 
water and chlorophyll 
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( 4 ) 


K = 0.040 + 0.0088(C/?/) + 0.054 (CM) 2 ' 3 
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Using equations (3) and (4), 
monthly MLD climatologic 
fields from NODC, and 
SeaWiFS Chl-a and PAR 
monthly products, we 
produced a monthly time 
series climatology of cross- 
shelf averaged (see Figure 1 
for transect location) Zm 
(MLD), Chl-a, and Zcr. The 
average was conducted from 
the coast to -340 km 
offshore to include the shelf 
and the MC upwelling zone. 
The seasonal time series is 
shown in Figure 7. 
Inspection of this figure 
reveals a classical behavior 
of the spring bloom. During 
winter, Zm exceeds or 
approaches Zcr, the time- 
integrated irradiance for 
photosynthesis is low, and 
despite elevated nutrients 
within the relatively deep 
MLD, production is low 
(Chl-a ~ 0.8 mg m' 3 on 
average). Around August- 
September, there is a fast 
decrease in Zm and a 
simultaneous large increase 
in Zcr, the mixed layer 
becomes shallower and still 
nutrient rich and the spring 
bloom occurs with peak 
production (Chl-a ~ 2.6 mg 
m"~ on average) around 
November. The mixed layer 
stays relatively shallow (< 
20 m) until February but the 
intense spring bloom consumes most of the nutrients by December-January, and a steep shallow 
thermocline prevents a steady nutrient input from below. Therefore, production quickly decreases 
and by autumn biomass returns to its background value of less than 1.0 mg m’ 3 of Chl-a. This 
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Figure 6. Monthly anomalies of SeaWiFS Chl-a (A, B, and C), AVHRR SST (D ? 
E, and F), and QuikScat wind speed and vectors (G, H, and I) for October, 
November and December 2004. QuikScat absolute monthly wind and vectors are 
shown at the bottom (J, K, and L). 
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scenario clearly shows that the Sverdrup critical depth model explains the mechanisms that drive the 
timing of the Patagonian shelf spring bloom. 



Figure 7. Seasonal time series of mixed layer depth (MLD. or 
Zm), Chl-a, and Sverdrup critical depth (Zcr). The thick grey line 
represents the average depth (-108 m) across the shelf and shelf 
break. 


Summer Autumn Winter Spring 3.4 Subsurface Cross-Shelf Transects 

Subsurface cross-shelf transects of 
monthly temperature, salinity, and 
temperature -derived nitrate for January 
(summer) and July (winter) 2004 are 
shown in Figure 8. The temperature and 
salinity fields originate from the HYCOM 
2003 - 2004 hindcast run, while the 
nitrate was obtained from the HYCOM 
temperature field and equation (1). 
Transects of monthly SeaWiFS Chl-a and 
MODIS SST are also shown. The red line 
on each of the subsurface fields is the 
monthly MLD, also obtained from the 
HYCOM run, while the blue line 
represents the Sverdrup critical depth. 
Note the large difference in temperature stratification between summer and winter. There is no 
vertical salinity stratification in winter or summer since there are no large freshwater sources along 
the Patagonian coast. The shelf is highly temperature-stratified during summer while the thermocline 
is very shallow at the PSB (~20 m based on the 11°C isotherm depth). Nitrate concentrations over 
the shelf range from 1 to 4 pM, with cross-shelf gradients indicating that the MC subsurface 
intrusion is the likely nitrate source. The MLD is less than 20 m across the whole transect. 
Conversely, during winter, the whole shelf is mixed almost all the way to the bottom, the 
thermocline over the shelf break is much deeper, and nitrate concentrations over the shelf range from 
12 to 20 pM. 

These HYCOM-derived cross-shelf parameters (T, S, N0 3 ) are reasonably similar to the 
values obtained by Carreto et al. [1995] from in situ data along a cross-shelf transect offshore of Mar 
del Plata (38°S - 40°S) for the same months, except for salinity being fresher nearshore in their data 
due to the proximity to the Rio de La Plata mouth. In addition, the derived nitrate using HYCOM 
temperatures does not include by definition the biological uptake during the spring bloom and thus 
overestimates summer nitrate within the mixed layer. The Carreto et al. [1995] transects revealed 
that a very strong and shallow (20 - 30 m) thermocline and nitracline developed during February 
and March 1989 with nitrate values much less than 1 pM above it. Chlorophyll subsurface maxima 
of 2.0 mg m’ 3 occurred at the nitracline depth. The chlorophyll stratification during the previous 
spring (October-November 1988) was quite different, with surface values across the shelf ranging 
from 2 to 7 mg m" 3 , and decreasing to less than 1 mg m’ 3 around 50 m. 

Note that the mixed layer depth ( Zm , red line) exceeds or is equal to the Sverdrup critical 
depth (Zcr, blue line) in July and Chl-a is relatively low across the shelf, while in January the Zcr far 
exceeds Zm and Chl-a is much higher (up to 6.0 mg m' 3 ), in accordance with the discussion given in 
the previous section. Also note that in January 2004 the minimum Chl-a occurs near 380 km 
offshore near the MC upwelling zone where the MLD is shallowest. This is probably a result of the 
faster nutrient depletion in surface waters by the spring bloom within the shallow MLD. However, it 
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has to be considered that satellite estimates correspond to chlorophyll in surface layers and, 
therefore, sub-surface peaks can occur which are not accounted for in our analyses. Note that, 
according to Figure 8, the nitrate isolines bend towards the surface near 380 Km offshore and this 
may have an effect on sub-surface phytoplankton growth. 



Figure 8. Left panel from top to bottom: January 2004 transects of SeaWiFS Chl-a and MODIS SST, 
and subsurface HYCOM temperature (7), salinity (5), and temperature-derived nitrate (N0 3 ) along 
the cross-section shown in Figure 1. The red line on the T, S. and NO ; cross-sections is the mixed 
layer thickness, while the blue line is the Sverdmp critical depth ( Zcr ). Right panel: same as left 
panel, except for July 2004. Note that the Chl-a and SST scales for January 2004 and July 2004 are 
different. 

To better illustrate the differences between the 2003 and 2004 bloom, Figure 9 shows the 
same fields, except for the addition of PIC, shown in Figure 8 but for December 2003 and 2004 
instead. The most striking difference between the two years is that during 2003 the bloom occupied 
the entire shelf with Chl-a concentrations ranging from 2 to more than 12 mg m’ 3 , while in 2004 the 
shelf Chl-a concentrations were mostly below 1 mg m°, except in the MC upwelling zone where 
Chl-a exceeded 2 mg m' 3 . The opposite occurred with PIC concentrations that were much less than 2 
mmol C m' 3 in 2003 but greater than 20 mmol C m’ 3 in 2004 at the PSB. One possible, and perhaps 
most probable, explanation for these differences is the control of productivity by the irradiance- 
mixing regime. The MLD during December 2003 was on average 7 in deeper than in December 
2004. This is a reflection of 2°C colder SSTs, on average, occurring during December 2003. During 
December 2003 the SST ranged from 9.6°C to 17.3°C with a mean of 13.6°C, while in December 
2004 they ranged from 12.0°C to 18.6°C with an average of 15.6°C. Note also that the Zcr is much 
deeper than Zm, but more so during December 2004, which is not shown because it exceeds 100 m. 
As explained before, a critical condition for the onset of the spring bloom is Zcr » Zm. However, 
the shallower the mixed layer the faster the nutrients are consumed and new nutrient supply is 
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greatly reduced. Simple integration of available total nitrate on the shelf revealed that the difference 
between 2003 and 2004 was less than 10%, so the higher productivity during December 2003 cannot 
be accounted for by changes in apparent nutrient supply from the MC at the PSB. Although the 
dynamics of nutrient uptake and supply cannot be inferred from the available data to ascertain the 
role of nutrients in the 2003 bloom, vertical mixing and advection seem to be the most probable 
physical explanation. Note that the PIC concentration, a proxy for coccolithophore presence, is 



Figure 9. Left panel from top to bottom: December 2003 transects of SeaWiFS Chl-a and calcite (PIC), 
MODIS SST, and subsurface HYCOM temperature (7), salinity (5), and temperature-derived nitrate (N0 3 ) 
along the cross-section shown in Figure 1. The red line on the T , S, and N0 3 cross-sections is the mixed 
layer thickness, while the blue line is the Sverdrup critical depth ( Zcr ). Right panel: same as left panel, 
except for December 2004. Note that the Chl-a scale for December 2003 and 2004 are different. Also note 
that the Zcr does not show on the December 2004 (right panel) cross sections as it is deeper than 100 m. 


highest at the PSB during 2004 where the MLD is shallowest. This is in agreement with our previous 
reasoning regarding the succession of phytoplankton groups, e.g., diatoms and dinoflagellates 
followed by coccolithophores. 
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3.5 Iron Sources and Possible Transport Mechanisms Conducive to Elevated Biomass 

The origin of macro-nutrients, especially nitrate, that support the intense blooms in the 
Patagonian shelf and shelf break can be traced to the MC upwelling and subsurface intrusions onto 
the shelf. In fact, modeling studies show that MC upwelling can be fairly intense in the Patagonian 
shelf-break [Matano and Palma, 2008], It is also conceivable, and ecologically required, that 
sufficient dissolved and bioavailable Fe originates from the MC upwelling zone but the actual Fe 
concentrations in water remain unknown due to the lack of sufficiently dense in situ observations to 
sample the relatively narrow MC. However, comprehensive field studies [Niencheski et ah, 2007; 
Windom et ah, 2006] along the Patos-Mirim lagoon system revealed a large, previously 
unrecognized source of dissolved iron to the South Atlantic Ocean via submarine ground water 
discharge (30.5°S to 32.5°S). These studies revealed that the iron flux from these lagoons to the 
coastal ocean is equal to about 10% of the soluble atmospheric Fe flux to the entire South Atlantic 
Ocean, implying that its effect on the Patagonian shelf could be potentially large. However, how this 
source correlates to the 2003 anomalous bloom situation is difficult to determine. 

Shelf sediment resuspension may also contribute to dissolved Fe inputs to the water column 
as shelf sediments tend to accumulate organic and inorganic compounds that sink down from the 
euphotic zone. Also, the MC may entrain dissolved Fe near the shelf break east of the Malvinas 
Islands (Falklands), a region that has been previously identified as having relatively high dissolved 
Fe [Bowie et ah, 2002]. Note in Figures 1 and 4 that the bloom around the Malvinas Islands in 2004 
was much more intense in size and magnitude than in 2003. If the MC is an important Fe source, the 
bloom around the Malvinas in 2004 may have depleted the local dissolved Fe concentrations, thus 
significantly reducing its transport northward by the MC. 

Interannual changes in mixed layer depth, especially during spring-summer, MC upwelling 
and location of the BMC can all contribute to changes in the Patagonian shelf production and 
biomass. Another possible Fe source is the coastal region off the southern end of South America 
where the astronomical tides are very strong (promoting intense shelf mixing), and the rivers and 
glaciers discharge terrigenous dissolved materials onto the shelf. The strong MC advection may 
capture and transport these materials northward. 
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Figure 10a. SeaWiFS monthly 
composites of Chl-a for December 2003 
(left) and December 2004 (right). 
Equivalent HYCOM currents averaged 
over 0 - 40 m are superimposed. Note 
that the vector scales are different for 
depths (h) less and greater than 300 m. 
The black contour on the shelf between 
36°S and 44°S is the 5 mg m' 3 Chl-a 
contour for December 2003, which is 
repeated for December 2004 as a 
reference. The 300 m bathymetry contour 
is also shown. 
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Figure 10 compares the circulation and water masses present on the PSB and inner shelf 
between December 2003 and 2004. Figure 10a shows the monthly SeaWiFS Chl-a and upper 40 
meters HYCOM mean currents for the two Decembers. During December 2004, a deflection of the 
MC is clearly shown in the inner shelf from 37°S to 43°S. A MC branch carrying slope waters 
deflected inshore from the PSB and a region of Chl-a above 5 mg m' 3 (black contour) resided on the 
shelf between the shoreward limit of the deflection and the PSB. Figure 10a clearly shows that there 
was no MC deflection during December 2004 and the Chl-a concentrations within the same region 
of the shelf were below 0.5 mg m°. Although there is indirect evidence that these MC deflections 
may be topographically induced [Franco et al., 2007, in press], it is possible that the 2003 deflection 
was one of the strongest during the SeaWiFS period. However, since the HYCOM simulation is only 
available for the period of November 2003 - December 2004, it is not possible to ascertain the 
anomalous nature of the December 2003 MC deflection. 

Further evidence of the effect of the MC deflection on the shelf bloom is provided in Figure 
10b, which shows a T-S diagram using values obtained from HYCOM monthly averaged 
temperature and salinity at all available depths within the bounds of the Chl-a 5 mg m" 3 contour. The 
T-S points are shown for December 2003 (circles) and December 2004 (crosses). The labels SASW 
and STSW indicate the approximate T-S characteristics of the Subantarctic Shelf water and 
Subtropical Shelf Water, respectively, as described by Piola et al. [2000], During December 2003, 
the shelf region where the bloom developed was highly influenced by the colder and fresher SASW, 
while during December 2004 the same shelf region was predominantly influenced by the warmer 
and saltier STSW. The colder SASW appears to originate from the outer shelf and the western edge 
of the MC (Figure 10a) and thus carrying higher nutrient concentration. Thus, it is quite possible that 
an inshore deflection was a major supplier of nutrients for the anomalous 2003 bloom. The 
northward displacement of the BMC, combined with SASW shelf intrusion caused by the deflection, 
have therefore a significant impact on the interannual variability of biomass on the Patagonian shelf. 

Another mechanism capable of elevating the Patagonian shelf biomass is the atmospheric 
deposition of Fe, which can originate from desert dust or volcanic ash. This fertilization from above 
is much harder to demonstrate since it requires sophisticated measurements of particle fallout at the 
ocean’s surface during hard to predict dust or volcanic ash events. Another difficult aspect is to 
determine the quantity of soluble Fe present during a given fallout event. Even so, previous studies 
provided evidence that dust and volcanic ash may play a significant role in the Fe fertilization of the 
southern ocean, albeit uncertainties regarding this process still abound. 

At southern South America, there is only scarce aerosol elemental composition data based on 
continuous atmospheric monitoring. Once these data are more abundant, they can be used to identify 
associations between the semi-arid Patagonian desert dust emission and the Patagonian shelf-break 
primary production. Presently, only models and remote sensing data have provided evidence of a 
possible association between the two above parameters [Erickson et al, 2003; Gasso and Stein , 
2007], Dust measurements conducted in Patagonia, reported by Gaiero et al. [2003], have clearly 
shown a significant correlation between the total aeolian dust flux and the iron suspended in the local 
atmosphere (R =0.95). Nevertheless, the total flux of dust onto the South Atlantic Ocean remains 
uncertain. Modeled data have estimated a range from -0.20 mmol Fe/m /yr [Dace and Tindale, 
1991] to -10 mmol Fe/m 2 /yr [Tegen and Fung, 1994], Considerations on the dissolved fraction of 
iron (DFe), which is effectively bioavailable to the marine biota, are also based on scarce data: Duce 
and Tindale [1991] estimated that for the South Atlantic, DFe flux into the ocean would be -0.08 
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10 g/yr, considering a conservative model where 10% of the iron dissolves in the sea. In addition, it 
is believed that chemical transformation and solubilization of the iron may occur during the aerosol 
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form, along the transport of dust from the desert site to the ocean, due to photochemical reactions 
and/or the effect of the acidic atmospheric environment. In the remote Southern South American 
continent, over Patagonia, acidic effect is unlikely to occur regarding the long distances with respect 
to main South American industrialized regions. 

Nevertheless, the oxidized sulfur (S0 X =S02 + 
sulfate) has been deemed recently as one of the most 
important chemical atmospheric species in the 
Southern South America region due to its 
acidification properties and potential Fe solubility 
effects. The source stems from combined effects of 
the numerous active volcanoes along the Chilean 
Andes, the copper concentrate smelters, and the 
emissions of at least two megacities: Santiago with ~ 
7 million habitants and Buenos Aires with ~12 
million habitants. Unlike North America and Europe 
that started reduction of SO 2 emissions from fossil 
fuel combustibles since the 1970s, petroleum 
byproducts employed in South America are still 
enriched with sulfur. Just before the October - 
December 2003 bloom, during the austral winter 
months of August-September 2003, a small eruption 
with Volcanic Explosivity Index or VEI = 0-1, took 
place at the Nevados de Chilian Volcanic Complex 
(36°50'S) [Naranjo and Lara , 2004], This event was 
potentially able to compose a combined contribution 
of desert dust and volcanic ash clouds. Although the 
VEI value does not represent an intense volcanic 
event (VEI scale varies from 0 to 8 and correspond to 
plume height restricted to the boundary layer), two 
aspects may be taken into account: (1) the eastern 
predominance of winds (from Andes to the Atlantic 
Ocean) at the low levels, which is favorable to dust 
delivery into the South Atlantic; and (2) the relative 
proximity of the Andean Volcanic arc with the 
Patagonian Shelf-break. Although the event of 29 
August 2003 was of relatively low magnitude, 
volcanic plumes can travel short continental distances 
before reaching shore areas in the Atlantic Ocean. 
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Figure 10b. T-S diagram using temperature and 
salinity values from HYCOM within 0-50 m and 
inside the region bound by the 5 mg m' 3 shown m 
Figure 10a. Note that the nitrate (N0 3 ) 
concentration derived from Equation (1) is 
represented by the color scale. The labels SASW 
and STSW indicate the approximate T-S 
characteristics of the Subantarctic Shelf water and 
Subtropical Shelf Water, respectively, as described 
by Piola et al [2000]. 

Figure 11 shows daily dispersion of one-hour, near-tropospheric integration of total ash and plume 
trajectory inferred from NOAA’s HYSPLIT Model between 30 August 30 and 2 September 2003. 
The impact over the South Atlantic is evident, even though there is a lack of precise data on 
emission factors. 
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Figure 11. Plume dynamics inferred from NOAA HYSPLIT Model from 2003 eruption at Nevados de Chilian. Plume 
positions in A, B, C and D correspond to days: August 30 th , August 31 st , September 1 st and September 2 nd , respectively. 
Concentration scale represents a relative value based on total emission. 

Iron in volcanic ash (differently from the desert dust) is composed of both iron oxides Fe 2 C >3 
and FeO, the second form constitute the soluble Fe(II) and in many geochemical analysis of Chilean 
volcanic ashes, its relative abundance may reach 2-3 times higher than the first form, which is the 
insoluble form Fe(III). Also, volcanoes are sources of SO 2 , which can promote pH acidification of 
clouds and lead to desert iron solubilization, independently of a soluble iron volcanic source. 
Therefore, despite the lack of in situ verification, Fe enhancement via atmospheric deposition during 
the 2003 anomalous bloom cannot be ruled out. The anomalous nature of the bloom may have been a 
result of a combination of physical factors, including the deflection of the MC onto the shelf and Fe 
atmospheric deposition. 

3.6 Further Investigation of Physical-Biogeochemical Interactions using Analytical and 
Numerical Models 

Additional analyses consisted of: 

1) Use of higher time resolution datasets, 8-day for satellite data and daily for HYCOM. Also, 
products from the biogeochemistry/ecosystem/circulation (BEC) model of Moore et al. (2004) were 
used. 
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2) Evaluation of nutrient transport onto the shelf based on daily time series of HYCOM u, v velocity 
components and derived nutrients from daily HYCOM temperature ( T) and salinity (S, for silicate 
only). 

3) Time series of phytoplankton growth were obtained using two methods, satellite-derived and 
model-derived. 

3.6.1 Nitrate Onshore Flux 

Figures 12 and 13 show MODIS Chl-a, HYCOM currents, and HYCOM-derived nutrient 
distributions for November 2003 and September 2005, which correspond to the highest and lowest 
biomass along the cross-shelf transect during the peak of the annual phytoplankton bloom, 
respectively. Figure 14 shows the nutrient regression equations. Only the nitrate onshore flux was 
computed. 

In finite difference form, the components of the nutrient flux (F x , F y ) at a given model grid 
point are given by 


F x 

F y 


..... .a N 0 3 (/' + 1, j) - N0 3 (/ - 1, j) 
~ { ’ J) 2 Ax 

= r -s no 3 (/, y + 1 ) - no 3 (?, y - 1) 

W,J) 2 Ax 


( 5 ) 


where u and v are the velocity components, NO3 the nitrate concentration, dx and dy are the grid cell 
widths in the east and north directions. The onshore nitrate flux was computed by rotating the flux 
components 45°, the local orientation of the 100 meter isobath, to obtain the cross-isobath flux. 

Figure 15 shows a map of the study region indicating the along-shelf transect near the 100 
meter isobath where the onshore nitrate flux was calculated. The estimated annual mean onshore 
nitrate flux is 1.5x10" mol m' yr . Figure 16 shows time series of SeaWiFS Chl-a averaged along 
the cross-shelf transect, onshore nitrate flux, and mean mixed layer nitrate along the same transect 
(blue line in Figure 15). Inspection of Figure 16 does not reveal any significant increase in onshore 
nitrate flux during the peak of the 2003 anomalous bloom. 

3.6.2 Phytoplankton Carbon Biomass and Growth using the Carbon-based Production Model 

The following analysis is based on the Carbon-based Productivity Model (CbPM) described 
in Behrenfeld et al. [2005] and Westberry et al. [2008]. 

The satellite-derived phytoplankton carbon biomass (Csat) is given by 

C SAT (mgC nT 3 ) = SF [b bp (443) - b bp (443) NAP ] (6) 

The scaling factor ( SF = 13000 mg-C m' 2 ) was derived such that the resulting C is -25-40% 
of total POC as estimated using the same backscatter values. The value for &6 p(443)nap, 0.00035 m' 1 , 
which is the correction for contributions due to a background of nonalgal backscattering particles. 
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was estimated by type II regression of monthly mean satellite-derived values of 6^(443) and Chi -a 
(see Behrenfeld et al., [2005]). 

The satellite-derived phytoplankton growth (p sat ) is given by 

) = / / max (Chl/QX ~ \}~ e ? ] 

The constant ju m ax , assumed to be 2 d’ 1 , approximates the observed maximum growth rate for a 
natural phytoplankton population ( Banse , 1991). The physiological response to light 

(photoacclimation) is characterized by a decreasing exponential function of light. More specifically, 
it results from a least squares fit of the satellite-derived Chl.C distribution as a function of monthly 
median mixed layer light intensities, which results in the expression 

(Chi / C) N . Tmax = [0.022 + (0.045 -0.022)e" 3 *] (8) 

The above expression is a representation of Chl.C in nutrient-replete, optimal growth 
conditions, e.g., the maximum potential Chi: C for a given irradiance. The two constants in the above 
equation reflect a low-light maximum in Chl.C (0.045 mg ChU mg C) and a high-light asymptotic 
value (0.022 mg ChU mg C). 

Ig (model-photons m ' 2 d ' 1 IT 1 ) is the median mixed layer light level, obtained using SeaWiFS 
PAR and A 490 , and the spectral fractionation numerical method [Austin and Petzold, 1986; Voss, 
1992] 

where Ez m id is the spectral median mixed layer light level in moles-photons m ' 2 d ’ 1 nm ’ 1 given by 

E Zmld {A) = 0.915E o {A)e ( - A) ^ (10) 

2 1 1 

The surface spectral PAR E 0 (X), in moles-photons m" d' nm’ , is obtained from the 
following equations that fractionate the integrated surface PAR into nine band components using the 
spectral coefficients/^) in nm ’ 1 

E 0 (X) = fW PAR 

f(X) = [0.0029, 0.0032, 0.0035, 0.0037, 0.0037, 0.0036, 0.0032, 0.0030, 0.0024] (11) 

/i(nm)=[400, 412, 443, 490, 510, 555, 625, 670, 700] 
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Figure 12. Eight-day composite of MODIS Chl-a for November 9-16, 2003 (a), and HYCOM-derived 
nitrate (b), silicate (c), and phosphate (d) for November 15, 2003. HYCOM currents (white vectors) for 
November 15, 2003 are plotted over each image. 
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Figure 13. Eight-day composite of MODI S Chl-a for September 22-29, 2005 (a), and HYCOM-derived 
nitrate (b), silicate (c), and phosphate (d) for September 30, 2005. HYCOM currents (white vectors) for 
September 30, 2005 are plotted over each image. 
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Figure 14. (a) Regression of nitrate versus 
temperature: log 10 [NO 3 ]=-0.003999*T 2 + 0.009936*T 
+ 1.488. (b) Regression of phosphate versus 
temperature: log 10 [P0 4 ]=-0.0002051*T 3 + 
0.003571*T 2 - 0.05464*T + 0.4319. (c) Regression of 
silicate versus T/S: 

logio[Si03]=-16.38*(T/S) 3 + 17.91*(T/S) 2 -8.252*T 
+ 1.916, ^=0.87. 


The spectral diffuse attenuation, K ( /(l) in m' 1 , 
as a function of K 49 0 and the spectral diffuse 
attenuation for water Kd w , is given by 

^ (490) ] + (2)<12) 

M(X) and Kd w (2) are tabulated constants, 
interpolated to the 9 discrete wavelengths from the 
multi-spectral values given in Austin and Petzold 
[1986]. The day length dl (in hours) is obtained from 
a function of year-day and latitude. The input 
parameters are for the period of November 7, 2003 - 
September 7, 2007 at 8-day resolution and averaged 
within the shelf domain delimited by the 5 mg/m' 
Chl-a contour obtained from the SeaWiFS Chl-a 
December 2003 composite. The input data originate 
from SeaWiFS GSM Chl-a (mgCW/m 3 ) and hb p (in 
m' 1 ) 8-day composites, SeaWiFS PAR (moles- 
photons/m 2 /d) and K 490 (m" 1 ) 8-day composites, and 
HYCOM daily temperature and salinity to derived 
nutrient concentrations, and HYCOM mixed layer 
thickness, together with PAR and K 490 , to calculate 
Ezmld- 

Figure 17 shows the time series of GSM Chl- 
a, HYCOM mixed layer depth, and satellite-derived 
growth, carbon biomass, and net primary production 
(NPP) averaged within the shelf domain defined 
above. Although there is some interannual 
variability in the growth, the variability does not 
correlate with the large carbon biomass observed 
during the anomalous bloom. If we assume that the 
satellite-derived growth is accurate and captures the 
underlying changes in nutrient availability, then the 
explanation for the large biomass observed in 2003 
must be elsewhere (grazer pressure, different 
proportion of functional groups, other physiological 
explanation?). One limitation of this method is that 
there is no distinction of functional groups and no 
explicit assessment of iron limitation effects, which 
we attempt to address in the next section. 
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3.6.3 Growth of Diatoms and Small 
Phytoplankton Based on Analytical Approach 
and Numerical Model Results 

Using Chl-a, SST, and nutrients from model 
simulations (Chl-a, SST, NO 3 , PO 4 , Si 03 , 
dissolved iron DFe, and NH 4 ), together with 
analytical equations relating nutrient concentration 
and euphotic zone light penetration to growth, an 
estimate of growth for diatoms and small 
phytoplankton was obtained for the study region. 
Monthly climatologic surface nutrients, 
chlorophyll and temperature originate from the 
global ocean Biogeochemical Elemental Cycling 
(BEC) ocean model of Moore and Braucher 

[2008], This is a global three-dimensional marine 
ecosystem model with several key phytoplankton 
functional groups, multiple limiting nutrients, 
explicit iron cycling, and a mineral ballast/organic 
matter parameterization. The model is run within a 
global circulation model. The model grid (102x116 
grid points) has a longitudinal resolution of 3.6° 
longitude, and a latitudinal resolution which varies 
from 1 ° to 2 °, with higher resolution near the 
equator. The caveat is that in the southwest 

Atlantic the model does not perform really well. The bloom in that region tends to be too big, and 
the coarse resolution doesn't capture the Malvinas Current well. Therefore, since the model does not 
resolve finer scale features, such as the Malvinas upwelling zone along the shelf break, an 
interpolation scheme was developed which uses the coarser resolution model products within the 
upper 40 meters to derive nutrient distributions based on higher resolution satellite images of Chl-a 

and SST. Figures 18 through 24 show monthly maps of Chl-a, Fe, NO 3 , SiC> 3 , PO 4 , and NH 4 , 

respectively, subset for the region of interest. These model results were used to derive nutrient 
regressions. The results are shown in equations 13 and 14 and the coefficients are given in Table 1. 

logic ( x ) = c o + C J + c 2 7 ’ 2 + c 3 lo §io {Chl) + c 4 T log 10 (Chl) + c 5 [log 10 ( Chi )] 2 

log 10 (Fe) = d 0 +d l Mo + d 1 Mo 2 + d } Lon + d 4 Lat + d 5 T + d 6 T 2 +d 1 log 10 (67?/) + 
c/ 8 [log 10 (Oi /)] 2 

Figures 25 and 26 show the 8 -day composite images of MODIS Aqua Chl-a, and the derived 
nutrients (DFe, NO 3 , and Si 03 ) from satellite SST and Chl-a for two time periods, December 1 1-18, 
2003, and September 22-29, 2005. These two periods highlight extreme high and extremely low 
biomass on the Patagonian shelf during spring. Note in these figures that the southern portion of the 
region, where the MC originates, the macro nutrients concentrations are relatively elevated while 
DFe concentrations are the lowest in the region. The shelf region, however, shows higher DFe and 


(13) 

(14) 



Figure 15. Map showing the along-shelf transect (red 
line) at which the on-shelf nitrate flux (large black 
arrow) time series was calculated. The annual averaged 
nitrate flux was 1.5x10 s mol NO 3 m ' 2 yr' 1 . The blue 
line indicates the cross-shelf transect through the 
middle of the shelf bloom. 
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macro nutrient concentrations, while the subtropical gyre region is poor in macro nutrients but 
relatively rich in DFe. This scenario is consistent with the satellite-derived biomass distribution. 





2003 2004 2005 2006 2007 2008 2009 

Figure 16. Time series of (a) mean SeaWiFS Chl-a along the cross-shelf transect 
(blue line in Figure 1), (b) onshore component of nitrate flux, and (c) mean mixed 
layer nitrate along the cross-shelf transect. The grey line in (a) is the 8-day Chl-a , 
while the red line is the low-pass filtered data. The same applies for (b) and (c), 
except that the grey lines represent daily data. The vertical lines are markers for the 
peak Chl-a concentration. 
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The growth equations for diatoms (DA) and small phytoplankton (SP) can be written as [Signorini 
and McClain, 2003]: 


[iDA=F phy (E) min [N H „ .PO„„„,SiO, 1 „, Fe lim ] 
PSP-F^XE) min[N lta ,P0 4mi ,Fe J 


Table 1. Coefficient values for nutrient polynomial regressions (equations 13 and 14) based on BEC ocean model. 


X 

Co 

Cl 

c 2 

c 3 

c 4 

c 5 




2 

r 

no 3 

2.056 

-0.175 

-0.863 

0.00435 

0.101 

-0.530 

- 

- 

- 

0.73 

Si0 3 

2.127 

-0.09008 

0.198 

-0.00176 

-0.0353 

0.03759 

- 

- 

- 

0.84 

po 4 

-0.155 

0.09624 

-0.886 

-0.0086 

0.08475 

-0.350 

- 

- 

- 

0.86 

nh 4 

-0.604 

0.03512 

-0.470 

-0.0063 

-0.02017 

-0.639 

- 

- 

- 

0.65 


^0 

di 

d 2 

d 2 

d 4 

d s 

d 6 

d 7 

ds 

2 

r 

Fe 

-1.369 

0.222 

-0.01389 

-0.01936 

0.02904 

0.08857 

-0.00426 

-0.201 

-0.01051 

0.52 


The light-dependent function isF phy (E) = ju^aE /[ju^ +(aE) 2 J , where ju max =3.0 d' 1 is the 

2 i i 

maximum phytoplankton growth, a=0.25 (Wm‘ )’ d’ is the slope of the P-I curve, and E is the PAR 
at depth z in W m' 2 . 


The nutrient limitation functions for nitrogen, phosphate, silicate, and iron are 


Nlim NH 41im + N0 31im 


NH 41im = 


NFL 


N0 3Um = 


(NH,+X m ) 
NO, 


P^ 41im 




Fe = 

rc lim 


(NO, + K n ) 

PO. 

(po,+a: p ) 

SiP 3 

(Si0 3 +* s ) 
DFe 

(DFe+^ F ) 


[l-NH 41jm ] 


(16) 


where, K N , A^ N h, Kp, Ks, and Kp are the half- saturation constants for nitrate, ammonium, phosphate, 
silicate, and iron, respectively. Following Lancelot et al. [2000], we used the values K n =1.0 pM, Kp 
=0.2 pM, and Kp =1 .2 nM for southern ocean diatoms. Ks was set to 1 .0 pM for diatoms only and Kp 
=0.12 nM for small phytoplankton. Timmermans et al. [2004] conducted growth experiments with 
four large open ocean diatoms from the Southern Ocean in natural (low iron) Southern Ocean 
seawater with increasing Fe concentrations. Their results showed that growth rates increased three- 
to six-fold with increasing DFe, and the species with the smallest cells had the highest growth rates. 
They report Fe half-saturation constants (Kp) for growth ranging from 0. 19 to 1.14 nM. 
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The averaged euphotic zone PAR, E(Ze ), was analytically derived from the integral of the 
light attenuation exponential function 


E r 

E(Z eil ) = — 1-1 -e * parZ “ J , where Eo is the surface PAR 


7 K 

e?/ PAR 


(17) 


From Morel [1988], Z eu and the mean Xpar within Z e „ can be evaluated from Chi as 


: 38.0 ChE 


K VAR =0.l2lChl° 


(18) 



mixed layer depth, (c) satellite-derived growth using 
CbPM, (d) satellite-derived phytoplankton carbon using 
CbPM, (e) satellite-derived net primary production in the 
mixed layer (p sat x C sat ). 


Figure 27 shows time series of MODIS 
Aqua Chl-a and SST, derived DFe and macro 
nutrients (NO 3 , SiC> 3 , PO 4 , NH 4 ), and the 
analytical diatom and small phytoplankton 
growths. The vertical lines in Figure 17 mark the 
peak Chl-a concentrations on the series. Note that 
the diatom growth is always smaller (about Vi) 
than the small phytoplankton growth. The Chl-a 
concentration peaks near maximum Ez e and 
shallowest MLD. The Chl-a concentration begins 
to rise when DFe is at its maximum 
concentration, the MLD is deepest, and Ez e begins 
to increase. The maximum growth of diatoms and 
small phytoplankton occurs almost always during 
the maximum Chl-a concentration. Table 2 
tabulates Chl-a, SST, nutrient, and growth values 
for the periods during which Chl-a concentration 
reached a maximum value. The lowest 
concentration of all nutrients occurs around the 
middle of December 2003 when Chl-a attained 
the highest concentrations (~25 mg m' 3 ) of the 
entire time series. 

Eight-day images of spatial distribution of 
diatom and small phytoplankton growth are 
shown in Figure 28 for December 11-18, 2003, 
and September 22-29, 2005. HYCOM surface 
currents for the same periods are superposed and 


seem to provide good agreement with the spatial distribution of phytoplankton growth, and thus 
biomass. 
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Table 2. Biogeochemical parameter values at the occasion of the Chl-a peak (see Figure 27). Tabulated parameters are 
MODIS Chl-a (mg m’ 3 ), MODIS SST (°C), derived near-surface (Z^) dissolved iron (DFe in nM), derived surface nitrate 
(N0 3 in pM), silicate (Si0 3 in pM), phosphate (P0 4 pM)), ammonium (NH 4 pM), and diatoms (pDA) and small 
phytoplankton (pSP) growths from analytical formulation. The last column tabulates the averaged PAR within the 
euphotic layer. 


Date 

Chl-a 

SST 

DFe 

no 3 

Si0 3 

po 4 

nh 4 

jiDA 

USP 

Ep.\R 

1 1/20/2002 

7.1 

12.6 

0.47 

3.17 

4.41 

0.41 

0.020 

0.56 

1.29 

27.3 

12/14/2003 

24.5 

14.0 

0.27 

2.00 

1.91 

0.21 

0.001 

0.38 

1.00 

28.1 

9/25/2004 

5.0 

9.4 

0.67 

5.34 

14.64 

0.72 

0.090 

0.49 

1.01 

16.1 

9/25/2005 

3.2 

9.8 

0.75 

4.59 

12.87 

0.72 

0.102 

0.54 

1.07 

16.3 

11/12/2006 

5.7 

12.2 

0.52 

3.60 

5.70 

0.48 

0.035 

0.55 

1.25 

24.2 

12/14/2007 

6.7 

13.6 

0.39 

3.27 

3.24 

0.36 

0.021 

0.48 

1.23 

27.5 
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Figure 18. Seasonal maps of top-layer ( Z pd ) chlorophyll 
re-gridded at 0.5°x0.5° resolution from the 
Biogeochemical Elemental Cycling (BEC) ocean model 
[Moore and Braucher , 2008], 
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Figure 19. Seasonal maps of SeaWiFS chlorophyll. 
Compare with BEC chlorophyll maps in Figure 18. 
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Figure 20. Seasonal maps of top-layer (Z^) dissolved iron 
(DFe) re-gridded at 0.5°x0.5° resolution from the 
Biogeochemical Elemental Cycling (BEC) ocean model 
[Moore and Braucher , 2008]. 
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Figure 22. Seasonal maps of top-layer (Z pc j) silicate (Si0 3 ) 
re-gridded at 0.5°x0.5° resolution from the 
Biogeochemical Elemental Cycling (BEC) ocean model 
[Moore and Braucher , 2008], 
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Figure 21. Seasonal maps of top-layer (Z pc j) nitrate (N0 3 ) 
re-gridded at 0.5°x0.5° resolution from the 
Biogeochenncal Elemental Cycling (BEC) ocean model 
[Moore and Braucher , 2008]. 
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Figure 23. Seasonal maps of top-layer (Z pc j) phosphate 
(P0 4 ) re-gridded at 0.5°x0.5° resolution from the 
Biogeochemical Elemental Cycling (BEC) ocean model 
[Moore and Braucher , 2008], 
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Note from the growth maps that diatoms should be 
more abundant on the shelf and along the subtropical 
convergence where there are indications from this 
analysis that DFe concentrations are higher. Small 
phytoplankton growth is more widespread in the 
open ocean as they require less DFe to grow, but 
seems to be stronger along the shelf and shelf break 
in general. In addition, both DA and SP growths are 
much higher during December 2003 than September 
2005. 


3.6.4 High Resolution, Empirically Derived 
Seasonal Surface Nutrient Climatology 


Seasonal composites of MODIS Chl-a and 
SST are shown in Figure 29. Monthly composites of 
MODIS Chl-a and SST, combined with empirical 
equations 13 and 14, were used to calculate 9-km 
monthly surface nutrient distributions. The monthly 
surface nutrient fields were then binned to produce 
seasonal (quarterly) composites. Daily HYCOM 
surface vector currents, from the period of 
November 2003 - July 2008, were averaged into seasonal composites and superposed on all the 
seasonal maps. Seasonal maps of nitrate, dissolved iron, silicate, and phosphate are shown in Figures 
30 and 31. These empirical estimates are based on regressions using the coarse resolution BEC 
model [Moore and Braucher, 2008] simulations (see Figures 20 through 23) and therefore do not 
accurately capture nearshore sources of nutrients, such as the Rio de La Plata plume and the 
discharge from the Patos and Mirim lagoon system (~30°S to ~37°S along the coast). The Rio de La 
Plata plume is known to be a significant source of silicate [Ciotti et al., 1995], and there are field 
studies showing that the lagoons are a significant source of iron enriched submarine groundwater 
that seeps through sediment in the nearshore region [. Niencheski et al. , 2007; Windom et al., 2006]. 
However, the BEC ocean model includes sedimentary and mineral dust sources of dissolved iron. 
The sedimentary Fe source is incorporated as a constant flux of 2 pmol Fe m' d' from sediments at 
the bottom of the ocean grid whenever depth is less than 1100 m. Dust deposition is from the 
climatology of Luo et al. [2003], The model has recently been applied to quantify ocean 
biogeochemical sensitivity to variations in mineral dust deposition (iron inputs) [Moore et al, 2006], 
the feedbacks between denitrification and nitrogen fixation [Moore and Don ey, 2007], and the ocean 
biogeochemical response to atmospheric deposition of inorganic nitrogen [Krishnamurthy et al, 
2007], 

A close examination of Figures 30 and 3 1 reveals that the Malvinas Current transports high 
concentrations of nitrate, silicate, and phosphate northward along the Patagonian shelf break. 
However, these same southern waters are depleted in iron, as shown by the empirically derived 
seasonal maps of DFe in Figure 30 (right panel). This is consistent with the low Chl-a concentrations 
observed by ocean color satellites, as shown by the MODIS seasonal composites in Figure 29 (left 
panel). Dissolved iron concentrations in surface waters reach a maximum on the Patagonian shelf 
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BEC Ocean Model NH 4 (fiM) 

Figure 24. Seasonal maps of top-layer (Z pd ) 
ammonium (NH 4 ) re-gridded at 0.5°x0.5° resolution 
from the Biogeochemical Elemental Cycling (BEC) 
ocean model [Moore and Braucher. 2008], 
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waters in winter, as shown by the July-August-September composite in Figure 30, most likely a 
result of intensified water column iron sediment flux via vigorous mixing by storms and convective 
overturning. By summer (January-February-March composite in Figure 30), the surface iron is 
significantly depleted as a result of uptake during the spring bloom. 

4. Summary and Conclusions 

This study reveals that the Patagonian shelf and shelf break regions are characterized by 
significant interannual variability of phytoplankton productivity. The intensity and geographic extent 
of the phytoplankton blooms are well documented by ocean color satellites but the mechanisms that 
drive the large interannual variability are not entirely understood. The anomalous 2003 bloom 
discussed in this paper is a clear example of how large the year-to-year changes can be. These large 
changes in phytoplankton productivity must be associated with corresponding changes of nutrient 
supply but its sources are not entirely known. In particular, the soluble Fe sources are still elusive 
despite a large body of literature that identifies riverbome, airborne, and oceanic sources of 
bioavailable Fe controlling biological activity in this region. 

Regarding airborne sources, it is possible to identify potential atmospheric factors controlling 
iron delivery to the South Atlantic and possible association with Chl-a blooms: (1) the climatic 
regime (aridity level) of the Patagonian semi-arid desert that would impose differential dust 
emissions to the Southern Atlantic Ocean, both in yearly and in long-term time scales, as retrieved 
from ice-cores from Antarctica; (2) episodic volcanic contribution that may act in two ways: as a 
source of soluble iron and as an indirect factor, through the S0 2 supply that would promote 
solubilization of atmospheric iron blown from the desert surface; and, (3) integrated contribution of 
several industrialized sites of South America with population higher than 6 million habitants (Rio de 
Janeiro, Sao Paulo, Buenos Aires, Santiago and Lima), with sulfur emissions underestimated by 
current inventories. 

It is important to determine Fe fluxes and concentrations in this important biologically active 
region, as well as the overall balance of macro and micro nutrients that control the phytoplankton 
seasonal and interannual variability. The MC upwelling along the PSB has been previously 
identified as a major source of macronutrients but Fe has not been measured in the MC waters to 
unequivocally explain the PSB spring-summer bloom biogeochemical dynamics. Flowever, 
interannual changes in the summer MLD modulate the intensity and transition of functional groups 
on the shelf and PSB bloom regions and episodic MC intrusions (deflections) fertilize the shelf 
region with nutrient-rich waters. Also, the mechanisms that drive the timing of the PSB spring 
bloom can be explained using the Sverdrup critical depth theory. 

The core of the Malvinas Current (MC) is a conduit of macro-nutrient-rich waters originating 
in the southernmost region of the study domain. However, the MC core has low phytoplankton 
biomass compared with the Patagonian shelf and shelf break regions. Although dissolved iron (DFe) 
concentration observations are still scarce in this region, the abundance of macro-nutrients and low 
Chl-a concentrations in the MC core could be an indication of depleted DFe waters. Alternatively, 
much higher concentrations of Chl-a on the shelf seem to indicate the presence of waters with more 
abundant DFe. Model simulations [Moore and Braucher, 2008; Moore et al., 2002a; Moore et al., 
2002b; Moore et al, 2004] combined with a global compilation of DFe observations, predict 
column-integrated inputs of DFe iron fluxes from mineral dust (0.3-0. 6 mmol Fe/m /yr) and 
sediments (0. 6-1.0 mmol Fe/m /yr) in the Patagonian shelf region that are high relative to other 
global ocean regions. 
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Figure 25. MODIS Chl-a (a), and derived near-surface (Z p d) dissolved iron (b), nitrate (c), and silicate (d) for 
December 11-18, 2003. 
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Figure 26. MODIS C7?/-tf (a), and derived near-surface (Z p( j) dissolved iron (b), nitrate (c), and silicate (d) for 
September 22-29, 2005. 
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Figure 27. Time series of (a) MODIS Chl-a , (b) MODIS SST and HYCOM mixed layer depth (MLD), (c) derived 
dissolved iron (DFe) and mean euphotic depth PAR (E Ze ), (d) surface derived nitrate and silicate, (e) phosphate and 
ammonium, and (f) modeled diatom and small phytoplankton growth. All parameters were averaged within the 5 mg m' 3 
Chl-a contour obtained from the December 2003 composite on the shelf. The vertical lines are markers for the peak Chl- 
a concentration. 
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Figure 28. Near-surface (Z pd ), eight-day maps of diatom and small phytoplankton growth for 
December 11-18 2003 ( a and b) and September 22-29 2005. HYCOM surface currents for the same 
time period are superposed. 
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Figure 29. Seasonal composites of MODIS Chl-a (left panel) and SST (right panel). The white vectors represent 
seasonal HYCOM surface currents. 




Although this analysis provided further insights in the relationship between biomass and 
nutrient cycling, the driving mechanisms and DFe sources for the extremely elevated biomass on the 
Patagonian shelf region during the spring of 2003 remain elusive. Further investigations, including 
comprehensive field data collection and numerical modeling, with emphasis on iron sources and 
nutrient concentrations, will be indispensable to unravel the large interannual changes in biomass 
seen by ocean color satellites in the region. There are existing three dimensional, coupled physical- 
biogeochemical models for applications in continental shelf systems that can be configured for the 
Patagonian shelf region. One example is the Regional Ocean Modeling System (ROMS) coupled to 
a seven state variable (phytoplankton, zooplankton, nitrate, ammonium, small and large detritus, and 
chlorophyll) biological model, which has been used to study the nitrogen and carbon cycles of the 
US East Coast [Fennel et al., 2006; Hofmann etal., 2008]. 

The biological component of the model system is undergoing some upgrades, which include 
the inclusion of different sizes of phytoplankton and zooplankton, and dissolved organic matter 
(DOM). A 10-km resolution ROMS grid for the NorthEast North American (NENA) shelf was 
nested within the HYCOM North Atlantic (NA) basin grid. Nesting of the NENA domain within the 
NA model enables the capture of circulation features and variability in the regional domain that is 
forced in the NA basin. An analogous modeling approach could be used for the Patagonian shelf 
region. 
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Figure 30. Empirically derived seasonal composites of surface nitrate (N0 3 , left panel) and dissolved iron (DFe, right 
panel) using equations 13 and 14, and monthly MODIS Chl-a and SST. 
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Figure 31. Empirically derived seasonal composites of surface silicate (Si0 3 , 
panel) using equation 13, and monthly MODIS Chl-a and SST. 
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